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Q | Abstract 

, We consider the problem of thermoacoustic tomography (TAT), in which one needs to recon- 

struct the initial value of a solution of the wave equation from its value on an observation surface. 
We show that if some geometric rays for the equation do not intersect the observation surface, 
PLh ' the reconstruction in TAT is not Holder stable. 

< 

^ ■ 1 Introduction 

Thermoacoustic tomography (TAT) is a hybrid imaging method that combines the significant con- 
trast of electromagnetic and high resolution of ultrasound imaging. A biological object is irradiated 
£N) . by a brief electromagnetic (EM) pulse in visible light or radiofrequency range. Some fraction of the 

EM energy is absorbed by the tissues. Since EM absorption is much higher in tumors, knowing the 
distribution a(x) of the absorbed energy would provide some valuable diagnostic information. The 
■ absorption causes thermoelastic expansion in the tissues, and thus triggers a pressure (ultrasound) 

wave p(x,t) propagating through the body. The pressure is then measured on an observation surface 
S completely or partially surrounding the object. The initial pressure f(x) = p(x,0) is roughly pro- 
portional to the EM energy absorption a(x). One thus concentrates on the recovery of f(x) instead 
of a(x). 

The commonly accepted model of TAT is (e.g., [3, 24]) 



( p tt ( x , t) = c 2 (x) A p(x, t), x G H n , t > 0, 

{ p(x,0) = f(x), Pt (x,0)=0, (1) 

p(y,t) = g{y,t),y e s,t > o. 

Here c(x) is the the ultrasound speed at location x, g(y, t) is the pressure measured at location 
y € S and time t. We assume that the sound speed c{x) is smooth and there are Co,co > such 
that Co < c(x) < Co for all x € IR n . The main problem of TAT is: given c(x), reconstruct / 
from g. Significant progresses have been made in mathematics of TAT. The readers are referred to 
[1, 6, 12, 25] for thorough discussion and relevant references. 

In this text, we concentrate on the stability analysis of the reconstruction of / in a region of 
interest f2, which is a bounded open domain in IR n , from the measured data g on the observation 
surface S. This issue has been partially addressed in a number of papers [6, 10, 11, 14, 15, 17, 21, 26, 
27]. 

We first recall the analysis for the case of constant speed c(x) = c. It was shown that if the 
so-called visibility condition is satisfied, the reconstruction is only mildly unstable, similarly to the 
inversion of the standard Radon transform (see, e.g. [14], for more discussion and references). The 
visibility condition requires that for each point x in £1, every straight line passing through x intersects 
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the observation surface S. For instance, if S is a closed hypersurface surrounding f2 then the visibility 
condition is satisfied. 

Let us now consider the case of variable speed c(x). If f2 is enclosed by S, it was argued and 
demonstrated in numerics in [11] that if there are some geometric rays trapped inside f2, then the 
singularities of / that lead to these rays cannot be stably reconstructed. On the other hand, the 
authors of [21] recently proved that under visibility condition, which can be defined for variable speed 
using geometric rays instead of straight lines, the reconstruction is Lipschitz stable. In this paper, 
we prove a complementary result, which shows that the reconstruction is not Holder stable if the 
visibility condition does not hold. The visibility condition can be violated when either the data is 
incomplete (S does not completely surround 0, see, e.g., [14]), or the trapping phenomenon occurs 
(e.g., [11]). These two cases have the same instability, and we do not distinguish between them. 
Although we do not use this in the text, it should be mentioned that the Lipschitz instability can be 
obtained from the general framework proposed in [20]. 

The text is organized as follows. In section 2, we recall the notion of wavefront set, whose 
propagation is the central issue of stability analysis in TAT (as well as other types of tomography, 
e.g., [5, 8, 9, 15, 18]). We then present our instability result in section 3. 

2 Singularities and Wavefront set 

In tomography, singularities are usually related to sharp details, for instance the boundaries of 
objects, jumps in densities, or interfaces between tissues. In many cases, singularities (rather than the 
exact image) are of the main interest (e.g., they are the objects of reconstruction in local tomography, 
e.g. [13]). We will see that propagation of singularities also plays an important role in the stability 
analysis of reconstruction in TAT. 

We denote by V(\R n ) and V(\R n ) the standard spaces of test functions and distributions on IR". 
We now recall the definition of wavefront set (e.g., [23]): 

Definition 2.1 Let f G V'(M n ), x G M n , and £ € IR n \ {0}. Then f is microlocally smooth at 
(^Oj^o) if there is a function (p € Co°(IR n ) satisfying <p(xo) ^ and an open cone V containing £o> 
such that T(ipf) is rapidly decreasing in V. That is, for any N > there exists a constant Cn such 
that 

\Hvfm\<C N (0~ N , for allteV. 

Here T is the Fourier transform and (£) = (l + l^] 2 ) 1 ^ 2 . The wavefront set of f, denoted by WF(f), 
is the complement of the set of all (xq,^o) € IR™ x (JR n \ 0) where f is microlocally smooth. 

For example, if / is the characteristic function of an open set O with smooth boundary <9f2, then 
(xo,£o) ^ WF(f) if and only if xq € dQ and £o is perpendicular to the tangent plane T Xo dVL of d£l 
at xq. 

Propagation of the wavefront set of the solution p(x, t) for equation (1) can be precisely described 
(e.g., [16, 21-23]). Since the initial velocity pt(x, 0) is zero, each element (xq, £o) G WF(f) propagates 
in two opposite directions £o and — £o- Let us consider the bi-characteristics (x±(s), t±(s), £±(s), t±(s)), 
which are the solutions of the following Hamiltonian systems: 

±(s) = -c 2 (x)£, i(s) = T, 
£( s ) = ±Vc 2 (x)|£| 2 , f (s) = 
(x(0),t(0),e(0),r(0)) = (x o ,0,=f£o,c(z )|£o|)- 

Since c(x) is smooth and < cq < c(x) < Co, these bi-characteristics are well defined on s € 
IR+, and (x±(s),t±(s),£±(s),T±(s)) € WF(p) for all s. We denote the (x,t)-projections of these 
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bicharacteristics by 7£ + (xo,£o) an d TZ-(xq,^o). The following result can be easily obtained by basic 
tools of microlocal analysis (e.g., [23]): 

Theorem 2.2 Consider equation (1). Let O be an open subset of IR n x IR + . Assume that for all 
G WF(f), K + (x,£) andK-(x,C) do not intersect O. Then p G C°°{0). 

Let TZ x (x,^) be the x-projection of TZ+(x,£) U lZ-(x,£). TZ x (x,^) is, indeed, a connected smooth 
curve in IR n , which we call a (geometric) ray. The following result is a simple consequence of Theorem 
2.2: 

Corollary 2.3 Consider equation (1). Let V be an open subset o/IR n . Assume that for all (x,£) G 
WF(f), the rays K x (x,t) do not intersect V. Then p G C°°{V x IR+). 

3 Instability of reconstruction in thermoacoustic tomography 

Let us return to the main equation of TAT: 

ptt{x,t) = c 2 {x) Ap(x,t), VzGlR n , t>0, 

p(x,0)=f(x), pt(x,0)=0. U 

Let S be a closed piece of a hypersurface in IR n , T > 0, V = S x [0, T], and g be the restriction of p 
on T. We define the linear operator 

T:L 2 (Q) — » P'(T) 
/ — » 5- 

Here, we identify L 2 (0) with the subspace of L 2 (lR ra ) containing functions supported in Q. We 
assume that T is injective and consider the stability problem of the reconstruction of / from g. Let 
U and V be open sets in IR n such that U C ft and S CV. 

Theorem 3.1 Assume that there exists a nonzero vector £o € H n \ such that for all x £ U, the 
rays 7Z x (x,£,o) do not intersect V. Then the reconstruction of f from g is not Holder stable. That 
is, there do not exist \i > 0, 5 > 0, sq, si > 0, and C > such that 

Wfhnn) < C\\Tff HS0{r) , for all f G H s ^n) satisfying ||/|| H . 1(n) < 8. 

In order to prove this theorem, we need an auxiliary result. Without loss of generality, we can 
assume that G U and £o = (0, ..,0, 1). Since U is open, there is an open set Uo C IR n_1 and e > 
such that Uo x I C U, where / = (— e,s). Let us fix a nonzero function /o G Cq°(Uo). For each 
x G H n , we write x = (x',x n ), where x' G IR ra_1 and x n G H. We now consider 

X = {/ G L 2 (U) : f(x) = f (x')h(x n ), h G L 2 (IR), supp(h) C 7}. 

Then, X is an infinite dimensional closed subspace of L 2 (Q). 

Lemma 3.2 For all s > 0, T induces a linear bounded operator 

T\ X :X^ H S (T). 

In what follows, we first prove Lemma 3.2, and then show how to obtain Theorem 3.1 from it. 
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3.1 Proof of Lemma 3.2 



Let / G X, then WF(f) C U x {A£o : A / 0}. Since K x (x,\£ ) = K x (x,£ ), one has Tl x (x,£) = 
7Zx(x,t;o) for all G Assuming the condition in Theorem 3.1, one deduces that TZ x (x,^) 

do not intersect V for all (x,£) G WF(f). 

We consider equation (2) and let V(f) = p\v X R + - Due to Corollary 2.3, V\x is a linear operator 
from X to C°°{V x IR+). Let Vb be an open set in H™ such that S CV CV CV. Then, for all 
s > and T > 0, V induces a linear operator: 

V\ x ■ X^H S (V x [0,T]). 

We now prove that this operator is bounded. Indeed, we first show that V : L 2 (Q) — > L 2 (Vb x 
[0, T]) is bounded. Let u(.,t) = J^p(.,s)ds, then, from (2), u solves the equation 

f u tt = c 2 (x)Au(x,t), x G IR n , t > 0, , . 

\ u(x,0) =0, ut(x,0) = f(x). [6) 

Denoting E(t) = \\c~ 1 ut(-, *) Hz, 2 (iR Tl ) *) lli 2 (^)' one obtains E 1 is independent of t, due to 

the conservation of energy (e.g., [2]). That is, E(t) = E(0) for all t G IR+. Since u t (.,t) = p(-,t), 
u(.,0) = 0, and ut(.,0) = f, one derives 

l|c -1 p(-.*)lli2(R») + ll V <>*)ll]R™ = I|c _1 /IIl2(]r™)- 
Using the inequalities < Co < c(x) < Co, one deduces that there is a constant ^4 > satisfying 

\\p(;t)\\LHJR") < Mfh^, for all / G L 2 (0). 

Hence, 

lbll.L 2 (Vox[0,T]) < l|p(-^)llL 2 (IR n x[0,T]) < ^11 / llz, 2 (IR") > 

which proves the boundedness of V : L 2 (ft) L 2 (Vb x [0,T]). 

Since H s {Vq x [0,T]) is continuously embedded into L 2 {Vq x [0, T]), applying Propositions 3.4 
(see appendix), one concludes that V\x '■ X — > -fP(Vb x [0, T]) is bounded. 

We are now ready to prove Lemma 3.2. Since the restriction 7Z(p) = p\r, as a linear operator 
U : H S (V x [0, T]) — > H 8 ' 1 / 2 ^), is bounded for any s > \, we have 

T\x=KoV\x :X^H S {T) 

is bounded for any s > 0. 

3.2 Proof of Theorem 3.1 

We first recall some facts concerning singular value decomposition (e.g., [4, 7]). Let Hi, H2 be Hilbert 
spaces and A : Hi — > be a bounded injective operator. Let A* be the adjoint operator of A and 
B = A* A. Then B : i^i — > H2 is a positive definite bounded operator. Let us denote by {s 2 } the 
eigenvalues of B and by {e^} the corresponding unit norm eigenvectors. Then {ej}j is an orthogonal 
basis of Hi. Denoting fj = j-A(ej), it is simple to show {fj}j is an orthonormal set in H2. If 
range(A) is dense in H2, then {fj}j is an orthonormal basis of H2- 

For an operator A, we denote by {sj(A)}j the set of above sj, which are chosen to be positive 
and decreasing. Then {sj(A)}j are called s-numbers (or singular values) of A. We will need the 
following asymptotic behavior of s-numbers: 

Lemma 3.3 (e.g., [4, page 119]) 
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1. Let si > S2 and i\ : H S1 (F) — ► H S2 (T) be the natural embedding. There exists a constant 

s 2~ s l 

c\ > independent of j, such that Sj(i\) < c\j » . 

2. Let s > 0, e > 0, and %2 : Hq(—e,e) — > L 2 (—e,e) be the natural embedding. There exists a 
constant c 2 > independent of j, such that Sjfa) > Ci3~ s ■ 

Proof of Theorem 3.1 Suppose that there exist \i > 0, S > 0, so, si > 0, and C > such that 

Wfh^n) < C\\Tf\r HS0(r) 

for all / G H S1 (Q) satisfying ||/||jj s i(n) < Then for any / G H S1 (£1), applying this inequality to 
ii .M s f , we have 

L 2 (n) <c ( ll T /ll# s °(r) V 



lH s i(Sl) \ \\J \\H s i (SI) 

where Ci is a constant independent of /. That is, 



^(nj^CillT/II^^H/ll^n), (4) 

for all / G H Sl (Q). Let us consider the following subspace of space X from Lemma 3.2: 

X S1 = {/ : f(x) = f(x',x n )=u (x')h(x n ),he H S \I)}. 

We then have X Sl C X Pi H Sl (Q). We now prove that (4) does not hold true for all / G X Sl . Indeed, 
due to Lemma 3.2, for any / G X Sl C X, one has 

(5) 

Here, s > and C s is a general constant depending on s. Combing (4) and (5), we have 
r/ll£ s{ r)ll/lfef n) < C S ||/|| L2(Q) < C s ||r/||^ 0(r) ||/||^ (n) . 

That is, 

mn) Y~\ c J\Efh^Y . (6) 



Since /(x) = fo(x')h(x n ), direct computations show 



Inequality (6) gives 



/IIl 2 (S1) — 11/01^2(^-1)11^11^2(7), 
H«i(fi) < ll/o||jfn(]R™-i)||^||ifn(7). 



L2(/) < ^f '™r° (r) V^ fOT al1 h e Jtf>(l). (-) 



J "~ \ \\Tf\\Hs(T) 



Let s > so, we now prove that this equality implies [sjfe)] 1 M < C s [sj(ii)] M for all j, where zi,Z2 
are the natural embeddings h : H S (T) ^ iJ s °(T) and i 2 : ^(1) ^ L 2 (I). Indeed, let be the 

orthonormal basis of Hq 1 (I) such that 

IIMl 2 (/) = s j( i 2)\\hj\\ H s^ I y 
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Let a = («i, aj) € IR 7 and h = J2k<j a khk- Since {hj}j are also orthogonal in L 2 (I), 

W h Wh(I) = Y,\ a k\ 2 \\ h k\\l^I)=Y,\ a k\ 2 4(i2)\\h k \\ 2 H s 1(I) 

k<j k<j 

> s |fe)El^l 2 ll^ll^i(/) = ^fe)l!^ll^ 1( /)- 

That is, 

INIl2 (/) > Sj(i 2 )||fe||H»i(j). (8) 
Now, consider {<7j}j to be the orthonormal basis of H S (T) such that 

llffill^o(r) = Sj(ii)\\gj\\H°(T)- 

Fixing j, we can choose a = (a±, ...,otj) ^ such that, for h = J2 k <j a khk and f(x) = fo(x')h(x n ), 
T(/) belongs to the orthogonal complement of span{gi, ...,gj-i}. Since the embedding i\ : H S (T) — > 
H S °(F) has dense range, {gu}k is also an orthogonal basis of H S °(F). One then obtains Tf = 
Ek>jPkgk in both H S °(T) and H S (T). Therefore, 

r/nW) = E/**iiW> = Y,Pk4(h)\\g k \\h(r) < #i)r/nW)- 

That is, 

WTfWnsoiD^SjiiiWfWHspy (9) 
Combining (7), (8), and (9), one arrives at [sj(«2)] 1_M < Cs[sj(ii)] /1 , where C s is a constant inde- 

/ -. \ SQ — S 

pendent of j. Due to Lemma 3.3, this implies [j~ Sl } ( ^' < C s [j~ ~ ] M . Choosing s big enough such 
that si(l — fi) < — we have a contradiction by letting j — > oo. The proof is completed. I 

Appendix 

Proposition 3.4 Assume that X,Y,Z are Banach spaces and Z is continuously imbedded into Y. 
Let T : X — ► Y be a bounded linear operator such that T(X) C Z . Then T : X — ► Z is also 
bounded. 

Proof Due to the closed graph theorem [19], if it is sufficient to prove that T has closed graph. 
That is, if {xk}k C X such that x k — > x in X and Tx,t — > z in Z, then Tx = z. In fact, since the 
imbedding Z — > Y is continuous and Tx^ — > z in Z, we have Tx& -> z in 7. On the other hand, 
since x& — > x in X and T : X — ► V is continuous, we have Tx& — ► Tx in V. Due to the uniqueness 
of limit of the sequence {Txk} k in Y, we have Tx = z. This finishes the proof. I 
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